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We calculate the probability of recoilless emission and detection of neutrinos 
■ (Mossbauer effect with neutrinos) taking into account the boundedness of the par- 

ent and daughter nuclei in the neutrino source and detector as well as the leptonic 
mixing. We show that, in spite of their near monochromaticity, the recoillessly emit- 
ted and captured neutrinos oscillate. After a qualitative discussion of this issue, we 
£Sj | corroborate and extend our results by computing the combined rate of v e produc- 

tion, propagation and detection in the framework of quantum field theory, starting 
from first principles. This allows us to avoid making any a priori assumptions about 
the energy and momentum of the intermediate-state neutrino. Our calculation per- 
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mits quantitative predictions of the transition rate in future experiments, and shows 
^ c*j that the decoherence and derealization factors, which could in principle suppress 

neutrino oscillations, are irrelevant under realistic experimental conditions. 
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Soon after the discovery of recoil-free emission and absorption of gamma rays by Mossbauer 
in 1958 [H, Hj], it has been suggested by Visscher that a similar effect should also exist for 
neutrinos emitted in electron capture processes from unstable nuclei embedded into a crystal 
lattice [3|. In the 1980's, the idea was further developed by Kells and Schiffer 0, 0], who 
showed that bound state beta decay [6| could provide an alternative recoilless production 
mechanism. In this case, an antineutrino with a very small energy uncertainty would be 
emitted, which could then be absorbed through induced orbital electron capture [7|. Recently, 
there has been a renewed interest in this idea, inspired by two works by Raghavan jl, @] , in 
which the feasibility of an experiment using the emission process 



and the detection process 



5 H -> 3 He + e~ (bound) + v e 



3 He + e" (bound) + v e -> 3 H (2) 
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has been studied. The 3 H and 3 He atoms were proposed to be embedded into metal crystals. 
The detection process would then have a resonance nature, leading to an enhancement of the 
detection cross section by up to a factor of 10 12 compared to the non-resonance capture of 
neutrinos of the same energy. If such an experiment were realized, it could carry out a very 
interesting physics program, including neutrino detection with 100 g scale (rather than ton or 
kiloton scale) detectors, searching for neutrino oscillations driven by the mixing angle 6*13 at a 
baseline of only 10 m, determining the neutrino mass hierarchy without using matter effects, 
searching for active-sterile neutrino oscillations and studying the gravitational redshift of 
neutrinos @, |, M, El 



In this paper we consider recoillessly emitted and captured neutrinos - which we will call 
Mossbauer neutrinos - from a theoretical point of view. In our discussion, we will mainly 
focus on the 3 H- 3 He system of Eqs. (OQ) and (j2J), but most of our results apply also to other 
emitters and absorbers of Mossbauer neutrinos. 

One of our main goals is to resolve the recent controversy about the question of whether 
Mossbauer neutrinos would oscillate. It has been argued 12] that the answer to this question 
depends on whether equal energies or equal momenta are assumed for different neutrino 
mass eigenstates - the assumptions often made in deriving the standard formula for the 
oscillation probability. Moreover, a possible inhibition of oscillations due to the time-energy 
uncertainty relation has been brought up To come to definitive conclusions regarding the 
oscillation phenomenology of Mossbauer neutrinos, we employ a quantum field theoretical 
(QFT) approach, in which neutrinos are treated as intermediate states in the combined 
production - propagation - detection process and no a priori assumptions on the energies or 
momenta of the different neutrino mass eigenstates are made. 

We begin in Sec. |2] by qualitatively discussing how the peculiar features of Mossbauer 
neutrinos, and in particular their very small energy uncertainty, affect the oscillation phe- 
nomenology. We argue that oscillations do occur, and that the coherence length is infinite if 
line broadening is neglected. We then proceed to quantitative arguments in Sec. [3] and discuss 
a formula for the v e survival probability in the quantum mechanical intermediate wave packet 



formalism [14j . in which the neutrino is described as a superposition of three wave packets, 
one for each mass eigenstate. In Sec. HI we derive our main result, the rate for the combined 
process of neutrino production, propagation and detection in the QFT external wave packet 
approach. In this framework, the neutrino is described by an internal line in a Feynman 
diagram, while its production and detection partners are described by wave packets. Also 
in this section, for the first time, we calculate the rates for beta decay with production of a 
bound-state electron and for the inverse process of stimulated electron capture in the case of 
nuclei bound to a crystal lattice. We distinguish between different neutrino line broadening 
mechanisms and concentrate on the oscillation phenomenology, paying special attention to 
the coherence and localization terms in the v e survival probability and to the Mossbauer 
resonance conditions arising in each case. In Sec. El we discuss the obtained results and draw 
our conclusions. 



2. MOSSBAUER NEUTRINOS DO OSCILLATE 



Mossbauer neutrinos have very special properties compared to those of neutrinos emitted 
and detected in conventional processes. In particular, they are almost monochromatic because 
they are produced in two-body decays of nuclei embedded in a crystal lattice and no phonon 
excitations of the host crystal accompany their production, which ensures the recoilless nature 
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of this process. Therefore the width of the neutrino line is only limited by the natural 
linewidth, which is the reciprocal of the mean lifetime of the emitter, and by solid-state 
effects, including electro mag netic interactions of the randomly oriented nuclear spins, lattice 
defects and impurities jl, H, 16|, 13]. For 3 H decay, the natural linewidth is 1.17-10 -24 eV, but 



it has been estimated that various broadening effects degrade this value to an experimentally 
achievable Mossbauer linewidth of 7 = O(10~ u eV) ltj ]. Compared to the neutrino 
energy in bound state 3 H decay, E = 18.6 keV, the achievable relative linewidth is therefore 
of order 10 -15 . 

In the standard derivations of the neutrino oscillation formula it is often assumed that the 
different neutrino mass eigenstates composing the produced flavor eigenstate have the same 
momentum (Ap = 0), while their kinetic energies differ by AE ~ Am 2 /2E. For bound state 
tritium beta decay (CQ) and Am 2 = Am| x ~ 2.5 x 1CT 3 eV 2 one has AE ~ 7 x 10~ 8 eV, which 
is much larger than 7. One may therefore wonder if the extremely small energy uncertainty 
of Mossbauer neutrinos would inhibit oscillations by destroying the coherence of the different 
mass eigenstates of which the produced v e is composed. Indeed, if neutrinos are emitted 
with no momentum uncertainty and their energy uncertainty (~ 7) is much smaller than the 
energy differences of the different mass eigenstates, in each decay event one would exactly 
know which mass eigenstate has been emitted. This would prevent a coherent emission of 
different mass eigenstates, thus destroying neutrino oscillations. If, on the contrary, one 
adopts the same energy assumption, the momenta of different mass eigenstates would differ 
by Ap ~ Am 2 /2p, which would not destroy their coherence provided that the momentum 
uncertainty of the emitted neutrino state is greater that Ap; in that case, oscillations are 
possible. 

It is well known that in reality neither same momentum nor same energy assumptions are 



correct [18|, [19|, [20|, [21|, [22|; however, for neutrinos from conventional sources both lead to 
the correct result, the reason being that neutrinos are ultra-relativistic and the spatial size 
of the corresponding wave packets is small compared to the oscillation length. 1 The above 
assumptions are thus just shortcuts which allow one to arrive at the correct result in an easy 
(though not rigorous) way. However, Mossbauer neutrinos represent a very peculiar case, 
which requires a special consideration. 

Let us discuss the issue of coherence of different mass eigenstates in more detail. If one 
knows the values of the neutrino energy E and momentum p with uncertainties &e and a p , 
from the energy-momentum relation of relativistic particles E 2 = p 2 + m 2 one can infer the 
value of the squared neutrino mass m 2 with the uncertainty 



(2Ea E ) 2 + (2pa p ) 2 , (3) 

where it is assumed that oe and a p are independent. By oe and a p we will now understand the 
intrinsic quantum mechanical uncertainties of the neutrino energy and momentum, beyond 
which these quantities cannot be measured in a given production or detection process; a m i is 
then the quantum mechanical uncertainty of the inferred neutrino squared mass. A generic 
requirement for coherent emission of different mass eigenstates is their indistinguishability: 
the uncertainty a m 2 has to be larger than the mass squared difference Am 2 23|. From the 



above discussion, we know that for Mossbauer neutrinos corresponding to the 3 H- 3 He system 
one has Eoe ~ 10~ 8 eV 2 , which is much smaller than Am 2 ~ 10~ 3 eV 2 . Thus, whether or 
not Mossbauer neutrinos oscillate depends on whether or not 2pa p > Am 2 . 



1 It is also essential that the energy and momentum uncertainties of these neutrinos are of the same order. 



While the energy of Mossbauer neutrinos is very precisely given by the production process 
itself, this is not the case for their momentum. The neutrino momentum can in principle be 
determined by measuring the recoil momentum of the crystal in which the emitter is embed- 
ded. The ultimate uncertainty a v of this measurement is related to the coordinate uncertainty 
a x of the emitting nucleus through the Heisenberg relation a v a x > 1/2. Therefore, for the 
momentum uncertainty to be small enough to destroy the coherence of different mass eigen- 
states, 2pa p < Am 2 , the coordinate uncertainty of the emitter must satisfy a x > 2p/Am 2 . 
This means that the emitter should be strongly de-localized with the coordinate uncertainty 
a x of order of the neutrino oscillation length L osc = 4np/Am 2 ~ 20 m. This is certainly 
not the case, because the coordinate uncertainty of the emitter cannot exceed the size of the 
source, i.e. a few cm. In fact, it is even much smaller, because in principle it is possible to 
find out which particular nucleus has undergone the Mossbauer transition by destroying the 
crystal and checking which 3 H atom has been transformed into 3 He. Thus, a x is of the order 
of interatomic distances, i.e. a v ~ 10 keV, so that 



This means that Mossbauer neutrinos will oscillate. The condition (J4]) is often called the 
localization condition, because it requires the neutrino source to be localized in a spatial 
region that is small compared to the neutrino oscillation length L osc . 

It should be noted that for the observability of neutrino oscillations the coherence of the 
emitted neutrino state is not by itself sufficient; in addition, this state must not lose its 
coherence until the neutrino is detected. A coherence loss could occur because of the wave 
packet separation. When a neutrino is produced as a flavour eigenstate, the wave packets of its 
mass eigenstate components fully overlap; however, since they propagate with different group 
velocities, after a time t coh or upon propagating a distance L coh ~ t coh , these wave packets 
separate to such an extent that they can no longer interfere in the detector, and oscillations 
become unobservable. The coherence length L coh depends on the energy uncertainty oe of 
the emitted neutrino state and becomes infinite in the limit oe — ► 0. 

From the above discussion it follows that the oscillation phenomenology of Mossbauer 
neutrinos should mainly depend on their momentum uncertainty, whereas their energy un- 
certainty, though crucial for the Mossbauer resonance condition, plays a relatively minor role 
for neutrino oscillations. Therefore, the equal energy assumption, though in general incor- 
rect, should be a good approximation when discussing oscillations of Mossbauer neutrinos. 
Adopting this approach, i.e. assuming the neutrino energy to be exactly fixed at a value E 
by the production process, one obtains for the i> e survival probability P ee at a distance L 



2pa p > Am 2 . 



(4) 





Here U is the leptonic mixing matrix, L° s £ are the partial oscillation lengths, 



AttE 



(6) 



L 



'jk - 




and the neutrinos are assumed to be ultra-relativistic or nearly mass-degenerate, so that 
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Eq. (jSj) is just the standard result for the v e survival probability. As expected, we do not 
obtain any decoherence factors if the neutrino energy is exactly fixed. We have also taken 
into account here that in real experiments the size of the source and detector are much 
smaller than the smallest of the oscillation lengths L° s £, so that the localization condition (Tj0) 
is satisfied. 



MOSSBAUER NEUTRINOS IN THE INTERMEDIATE WAVE PACKET 

FORMALISM 



Although Eq. (0) shows that neutrino oscillations are not inhibited by the energy con- 
straints implied by the Mossbauer effect, the assumption of an exactly fixed neutrino energy 
is certainly unrealistic. Therefore, we will now proceed to a more accurate treatment of 
Mossbauer neutrinos using an intermediate wave packet model [HI, 0, In this 

approach, the propagating neutrino is described by a superposition of mass eigenstates, each 
of which is in turn a wave packet with a finite momentum width. With the assumption of 
Gaussian wave packets, Giunti, Kim and Lee [HI, HI] obtain the following expression for the 
v e survival probability in the approximation of ultra-relativistic neutrinos: 



Pee = ^ \U ej \ 2 \U ek \ 2 exp 



2m- 



T osc 
L jk 



T coh 
L jk 



o^t- r OSC 

za p L jk 



Here 



2V2E 2 



- coh 

Jjk ~ <rp\*rn%\ 



(9) 



are the partial coherence lengths, a p being the effective momentum uncertainty of the neutrino 
state, and the oscillation lengths L°| c are given by Eq. E is the energy that a massless 
neutrino emitted in the same process would have, and the 0(1) parameter £ quantifies the 
deviation of the actual energies of massive neutrinos from this value. Since the energy 
uncertainty is very small for Mossbauer neutrinos, the mass eigenstates differ in momentum, 
but hardly in energy, so that £ should be negligibly small in our case. 

One can see that the first term in the exponent of Eq. (jHJ) is the standard oscillation phase. 
The second term yields a decoherence factor, which describes the suppression of oscillations 
due to the wave packet separation. For conventional neutrino experiments with non-negligible 
£, the third term implements a localization condition by suppressing oscillations if the spatial 
width a x = l/2a p of the neutrino wave packet is much larger than the oscillation length L°| c 
(cf. Eqs. (E]) and (pD)). However, we have seen that, due to the smallness of £, the intermediate 
wave packet formalism predicts this condition to be irrelevant for oscillations of Mossbauer 
neutrinos. 



4. MOSSBAUER NEUTRINOS IN THE EXTERNAL WAVE PACKET 

FORMALISM 

In the derivation of the quantum mechanical result discussed in the previous section, 
certain assumptions had to be made on the properties of the neutrino wave packets, in 
particular on the parameters a p and £. We will now proceed to the discussion of a QFT 
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Figure 1: Feynman diagram for neutrino emission and absorption in the 3 H- 3 He system. 



approach [26, 27, 28, 29|, 30, 31, 32, 33, 34, 35], in which these quantities will be automatically 
determined from the properties of the source and the detector. 

Our calculation will be based on the Feynman diagram shown in Fig. (TJ in which the 
neutrino is described as an internal line. We take the external particles to be confined by 
quantum mechanical harmonic oscillator potentials to reflect the fact that they are bound 
in a crystal lattice. Typical values for the harmonic oscillator frequencies are of the order of 
the Debye temperature 6^ ~ 600 K ~ 0.05 eV of the respective crystals (U 15]. Although 



this simplistic treatment neglects the detailed structure of the solid state lattice, it is known 
to correctly reproduce the main features of the conventional Mossbauer effect [36j , and since 
we are interested mainly in the oscillation physics and not in the exact overall process rate, 
it is sufficient for our purposes. As only recoil-free neutrino emission and absorption are of 
interest to us, we can neglect thermal excitations and consider the parent and daughter nuclei 
in the source and detector to be in the ground states of their respective harmonic oscillator 
potentials. 

In Sec. 14.11 we will develop our formalism and derive an expression for the rate of the com- 
bined process of Mossbauer neutrino emission, propagation and absorption. In Sees. 14.21 - 14.41 
we will then discuss in detail the effects of different line broadening mechanisms. 



4.1. The formalism 



Let us denote the harmonic oscillator frequencies for tritium and helium in the source 
by ^h,5 and ^Hc,s and those in the detector by cui^d, and Uu c ,d- In general, these are four 
different numbers because 3 H and 3 He have different chemical properties, and because their 
different abundances in the source and detector imply cl>h,s 7^ <^h,l> and uu c ,s 7^ ^h c ,d- We 
ignore possible anisotropics of the oscillator frequencies because their inclusion would merely 
lengthen our formulas without giving new insights into the oscillation phenomenology. The 
normalized wave functions of the ground states of the three-dimensional harmonic oscillators 
\^>a,b,o) are given by 



^A,B,o(x,t) 



m A UJA,B 



7i 



_, 3 

4 



cxp 



1 



-iE A ,Bt 



(10) 



where A = {H, He} distinguishes the two types of atoms and B = {S, D} distinguishes 
between quantities related to the source and to the detector. The masses of the tritium and 
3 He atoms are denoted by m H and mu e , and the coordinates of the lattice sites at which 
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the atoms are localized in the source and in the detector are xg and x^. The energies 
Ea,b of the external particles are not exactly fixed due to the line broadening mechanisms 
discussed in Sec. EJ but follow narrow distribution functions, which are centered around 
Ea,b,o — m A + \^a,b- For the differences of these mean energies of tritium and helium atoms 
in the source and detector we will use the notation 

Es,0 = -Eh,S,0 — -^Hc.S.O , Ed,0 = Eu,D,0 — -Efie.D.O • (11) 

Before proceeding to calculate the overall rate of the process of neutrino production, prop- 
agation and detection, we compute the expected rates of the Mossbauer neutrino production 
and detection treated as separate processes, ignoring neutrino oscillations. This calculation is 
very instructive, and we will use its result as a benchmark for comparison with our subsequent 
QFT calculations. 

The effective weak interaction Hamiltonians for the neutrino production and detection 
Hg and are given by Eqs. ( 1C2I) and ( 1C3I) of appendix C. We will first assume that 
the neutrino emitted in the recoil-free production process (CQ) is monochromatic, i.e. neglect 
the natural linewidth as well as all broadening effects. Likewise, we will neglect now the 
absorption line broadening effects in the recoilless detection process (T5]). A straightforward 
calculation gives for the rate of recoilless neutrino production 



T P = T X S , (12) 



where 



G 2 F COS 2 6 c 2 2 , 2 , „2 i 71 |2\ f Es,0 



\MR)\ 2 m 2 e {\M V \ 2 + g 2 A \M A \ 2 ) Kg (13) 



it \ m, 



2 



with Gf the Fermi constant, 8 C the Cabibbo angle, m e the electron mass, My and Ma the 
vector and axial- vector (or Fermi and Gamow- Teller) nuclear matrix elements and gA — 1-25 
the axial- vector coupling constant. Note that for the allowed beta transitions in the 3 H- 3 He 
system, My = 1 and Ma ~ v^3- The quantity ij) e (R) is the value of the anti-symmetrized 
atomic wave function of 3 He at the surface of the nucleus. The factor k$ takes into account 
that the spectator electron which is initially in the Is atomic state of 3 H ends up in the Is 
state of 3 He. It is given by the overlap integral of the corresponding atomic wave functions: 

us = J ^z=2,s( r T^z=i,s^)d 3 r 2 . (14) 

The factor X$ in Eq. ( TT21) is defined as 

1 X 



X s = 8 [rjs + — e V? = Y s e "*s , (15) 



where p = \/Eg — m? is the neutrino momentum 2 , and 



Vs= x , cr s = m K u; H ,5 + m Hc uj HCiS . (16) 



2 Since in this calculation we ignore neutrino oscillations, we also neglect the neutrino mass differences. 
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The energy spectrum p(E) of the emitted Mossbauer neutrinos in the considered approxima- 
tion is 

p(E) = r X s 8(E-E s ,o). (17) 
For the cross section of the recoilless detection process we obtain 

<r{E) = B X D 6(E-E Dfl ), (18) 

where 

B =4ttG 2 f cos 2 6 c \MR)\ 2 {\Mv\ 2 + 9 2 a\M a \ 2 ) k d . (19) 

The factor Kb here is defined similarly to in Eq. (fT4j) . Note that in the approximation of 
hydrogen-like atomic wave functions one has Kg = &d = 512/729 ~ 0.7. The factor in 
Eq. (1181) is defined similarly to the corresponding factor for the production process, i.e. 



_3 „2 



V _ V 

"3 



X D = 8 I n D + — ) e = y D e ^ D (20) 
Vd 



with 



Vd = x , cr D = m H w h ,d + m Ue u UCtD . (21) 

V m Hc ^He,D 

The Mossbauer neutrino production rate T p and detection cross section a(E) differ from those 
previously obtained for unbound parent and daughter nuclei respectively in Refs. j6| and 0] 
by the factors Xs and Xp. Note that in the limit mn u^s = ^He ^He,5, ^h,d = nine ^Ke,D, 
the pre-exponential factors Ys and Yd in Eqs. ( TT5l) and ( 1201) become equal to unity, so that 
Xs and Xp reduce to the exponentials, which are merely the recoil-free fractions in the 
production and detection processes (see the discussion below). 

For unpolarized tritium nuclei in the source the produced neutrino flux is isotropic; there- 
fore the spectral density of the neutrino flux at the detector located at a distance L from the 
source is p(E) / (AitL 2 ). The detection rate is thus 

p(E)a(E) dE = X S X D 8(E s ,o - E Dfi ) . (22) 



We see that it is infinite when the Mossbauer resonance condition E$fi = -Ed,o is exactly 
satisfied and zero otherwise, which is a consequence of our assumption of infinitely sharp 
emission and absorption lines. This assumption is certainly unphysical, and a realistic cal- 
culation should take into account the finite linewidth effects. We do that here by assuming 
Lorentzian energy distributions for the production and detection processes, which will be use- 
ful for comparison with the results of our subsequent QFT approach. In this approximation 
Eqs. ( fTTI) and ( |T8i) have to be replaced by 
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where 75- and 7^ are the energy widths associated with production and detection. The 
combined rate of the neutrino production, propagation and detection process is then 



4ttL 2 



p(E)a(E) dE 



To B 



XsX D 



(75 + 7d)/2tt 



4ttL2 (E s ,o - £d,o) 2 + (7s + 7d) 2 /4 

As can be seen from this formula, the Mossbauer resonance condition is 

(E s ,o - E Dfi f < (75 + 7d) 2 /4 • 



(24) 



(25) 



If it is satisfied, the neutrino detection cross section is enhanced by a factor of order 
(aZm e ) 3 '/\p e E e (js + Id)] compared to cross sections of non-resonant capture reactions 
u e + A —* A' + e + for neutrinos of the same energy (assuming the recoil-free fraction to 
be of order 1). For 75 + 7d ~ 1CT 11 eV the enhancement factor can be as large as 10 12 . 

We now turn to the QFT treatment of the overall neutrino production, propagation and 
detection process, first neglecting the line broadening effects. We derive the corresponding 
transition amplitude from the matrix elements of the weak currents in the standard way by 
employing the coordinate-space Feynman rules to the diagram in Fig. [TJ For the external 
tritium and helium nuclei, we use the bound state wave function ip A B,o( x , t) from Eq. (TTU]) . 
We obtain 



iA 



a Xi ati / a x 2 dt^ I I exp 



- -"lH^H,s|Xl 



-iEwsti 



( 



m Hc^He,S 
71 

^He^He,D 
71 

71 



exp 



exp 



-m Hc WHc,D|X2 - X D | 



-iE HeD t2 



exp 



3 



MW£\U, 



'ej I 



1 

2 

d A p 



m H o;H,D|X2 - x D 

ipo(fe-ti)+«p(x2-xi) 



(2tt)< 



rrij + le 



(1 + 7 5 )7„M eijD - 



(26) 



The Dirac spinors for the external particles are denoted by u a .b with A = {e, H, He} and 
B = {S, D}. Note that all spinors are non-relativistic, so that we can neglect their momentum 
dependence. The matrix elements Aig and Ai D encode the information on the bound state 
tritium beta decay and also on the inverse process, the induced orbital electron capture which 
takes place in the detector. They are given by 



G p cos 6 C 



MR) u Kc (M v 5% - g A M A <Ti 6?/V3)un n l/2 



S,D 



(27) 



The integrations over t\ and t 2 in Eq. (1261) yield energy-conserving 5-functions at the 
neutrino production and detection vertices. The spatial integrals are Gaussian and can be 
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evaluated after making the transformations Xi — > Xi + and x 2 — > x 2 + x^. We obtain 



d 4 p 



27i5(p - Es) 2n5(p Q - E D ) exp 



(2vr) 4 

M tx s M'S\U ej \ 2 u e , slfl {l-^) 



+ mj)e iph 



2 9 2 



- (1 + 7 & ) 7l ,u 



(28) 



where we have used the notation 

Eg = E^s — -^He,S j E D = Eft D — E^D , 

and introduced the baseline vector L = xo — xg. The quantity a p , which is given by 



1 



+ 



can be interpreted as an effective momentum uncertainty of the neutrino. Note that a, 



(29) 



(30) 



a 



P s 



a P D- We have also defined a constant 



J\f 



V 7T 



"lHe^He,,S 



7T 



7T 



( 



7T 



2tt 



2tt 



(31) 



containing the numerical factors from Eq. (1101) and coming from the integrals over Xi and 
x 2 . One of the 5-functions in Eq. (1281) can now be used to perform the integration over 
Po, thereby fixing p at the value p = Es = Ejj. To compute the remaining integral over 
the three-momentum p, we use a theorem by Grimus and Stockinger 30], which states the 
following: Let ^>(p) be a three times continuously differentiate function on M 3 , such that ip 
itself and all its first and second derivatives decrease at least as l/|p| 2 for |p| — > oo. Then, 
for any real number A > 0, 



d 3 p 



tp(p)e 



ipL 



|L|^oo 



iVAL 



0(L 



A - p 2 + ie L 
The validity conditions are fulfilled in our case, so that in leading order in 1/Lwe have 



(32) 



iA=-^MS(E s 



2L 



E 



cxp 



El 



M s M u j;\U ej \ 2 e l ^ E 's- m i L 



■ u e) slp,(l - 7 5 )(^- + mj)(l + 7 5 )7^ e ,D 



(33) 



where the 4-vector pj is defined as pj = (Es, (El 



2\l/2 



L/L). The Grimus-Stockinger 



theorem ensures that for L E^ 1 , where Eq is the characteristic neutrino energy, the 
intermediate-state neutrino is essentially on mass shell and its momentum points from the 
neutrino source to the detector. 

The transition probability V is obtained by summing \A\ 2 over the spins of the final 
states and averaging it over the initial-state spins. Note that no integration over final-state 
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momenta is necessary because we consider transitions into discrete states. The transition 
rate is obtained from V as V = dV/dT, where T is the total running time of the experiment. 
As we shall see, in the case of inhomogeneous line broadening V oc T for large T, so that 
T is independent of T in that limit. The same is true for the homogeneous line broadening, 
except for the special case of the natural line width, for which the dependence on T is more 
complicated (see Sec. 14.4ft . 



4.2. Inhomogeneous line broadening 

Inhomogeneous line broadening is due to stationary effects, such as impurities, lattice 
defects, variations in the lattice constant, etc. 1_5|, 113] - These effects are taken into ac- 
count by summing the probabilities of the process for all possible energies of the external 
particles, weighted with the corresponding probabilities of these energies. In other words, 
one has to fold the probability or total rate of the process with the energy distributions of 
tritium and helium atoms in the source and detector, pne,s(En Ct s), Ph,d(-Sh,d), Ph^-Eh^s) 
and pn e ,D( E Ke,D) ■ We obtain 



V 



dEft s dEftcS dEft e n dEu D 



■ PU,s(E HjS ) PHe,D(E}ie,D) PHe,s(En^ s ) P^d{E^d) |-4| 2 , 



(34) 



where |*4.| 2 is the squared modulus of the amplitude, averaged over initial spins and summed 
over final spins. Using the standard trace techniques to evaluate these spin sums and ne- 
glecting the momenta of the non-relativistic external particles, one finds 



V=T G% cos 4 6 c 



POO 

\MR)\ 4 El (\M v \ 2 + g 2 A \M A \ 2 ) 2 Y s Y D K S K D / dE u , s dE He , s dE He , D dE R , D 

Jo 



ttL 2 

6(E S — E D )p H: s(E U:S ) Puc,D(E UCtD ) pHe,s(EHe,s) Ph,d(-Eh,d) 



ej I 



\U ek \ 2 exp 



2E 2 



2 2 

m~ — mi 



2a 2 



(35) 



where Y$ and Yd were defined in Eqs. ( 1151) and (I2U1) . Here we have taken into account that 
for T 3> (Es — Ed)" 1 the squared (^-function appearing in l^] 2 can be rewritten as 3 



1 r T / 2 T 

[5(E S - E D )} 2 ~ — 5(E S - E D ) / dt e^~ E ^ = —5{E S - E D ) 

Alt J-T/2 ^ 



(36) 



3 The expression 5(Es — Ed) here should be understood as a 5-like function of very small width. For 
\E S - E D \ ~ 10~ n eV, the condition T > (E s - Ed)" 1 would require T > 1CT 4 s, which should be very 
well satisfied in any realistic experiment. 
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The overall process rate T is then obtained from Eq. (I35j) by simply dividing by T. Using 
the definitions of T and B given in Eqs. (|T5|) and (1T91 . one finds 



To B 



AttL 2 



YsY D / dE^s dEuc,s dE^ D dEu,D 



S(Es — E D )pn t s(Eii t s) Phc,d(Eho,d) PRe,s(En e ,s) Ph, d{Eu,d) 



|C/ efe | 2 exp 



2£^ - m) 



mi 



(37) 



Before proceeding to the computation of the remaining integrations over the energy dis- 
tributions of the external particles, let us discuss the expression in the last line in Eq. (I3T|) . 
In the approximation of ultra-relativistic (or nearly mass-degenerate) neutrinos, Eq. (J7J), 
the last exponential becomes the standard oscillation phase factor exp(— 2iriL/L° s ^) with 
the oscillation length defined in Eq. (El). The additional exponential suppression term 
exp[— (2 Eg — m| — mVi 1 2a 2 } is an analogue of the well-known Lamb-Mossbauer factor (or 
recoil-free fraction) 0, S, HJ, which describes the relative probability of recoil- free emission 
and absorption compared to the total emission and absorption probability. We see that for 
Mossbauer neutrinos this factor depends not only on their energy, but also on their masses. 



Therefore, if two mass eigenstates, 



and do not satisfy the relation |Am| fe | 



< (Tp, the 



emission and absorption of the lighter mass eigenstate will be suppressed compared to the 
emission and absorption of the heavier one. This can be viewed as a reduced mixing of the 
two states, which in turn leads to a suppression of oscillations. To stress this point directly 
in our formulas, we rewrite the corresponding factor as 



exp 



r (pT) 2 ' 






exp 


i . 


2a 2 

p 



(38 



where p™ n is the smaller of the two momenta of the mass eigenstates Vj and u/,, 



^min\2 



El 



max^m^ , 



m. 



(39) 



The first exponential in Eq. ( 1381) describes the suppression of the emission rate and the 
absorption cross section, i.e. is a generalized Lamb-Mossbauer factor, while the second one 
describes the suppression of oscillations. The condition |Am| fe | < 2o~ 2 enforced by this 
second exponential can also be interpreted as a localization condition: Defining the spatial 
localization <j x ~ l/2a p , we can reformulate it as L°| c > 4TTa x Es/<y p . Since the generalized 
Lamb-Mossbauer factor (the first factor in Eq. ( 1381) ) enforces Es < o~ p , this inequality is 
certainly fulfilled if \L° S £\ > 2na x holds. The latter, stronger^ localization condition is the 
one obtained in other external wave packet calculations (28^3 



34j and is also equivalent to 



the one obtained in the intermediate wave packet picture [24], ll4j and discussed in Sec. El 

Let us now consider the integrations over the spectra of initial and final states in Eq. (137]) . 
To evaluate these integrals, we need expressions for Pa,b, based on the physics of the inho- 
mogeneous line broadening mechanisms. To a very good app roximation, these effects cause a 
Lorentzian smearing of the energies of the external states 
are 



app 
371. 



37(, so that the energy distributions 



(40) 
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where, as before, A = {H, He}, B = {S, D} and Ea,b,o = m A + \^a,b- After evaluating the 
four energy integrals in Eq. (1371) (see appendix [X] for details), we obtain 



47T-L 2 27T 



l^efel 2 exp 



J'.* 
1 



(pT) 2 



exp 



2a 2 
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E 



5,0 



R„ _|_ 75 -ID 



+ 



5,0 



^d,o =F i 



7S+7D 



(41) 



In deriving this expression we have used the fact that the generalized Lamb-Mossbauer factor 
is almost constant over the resonance region and can thus be approximated by its value at 
E = 7;(E S) o + E D fl). The quantities in Eq. (14T|) are given by 



A jk 



exp 



Am 2 jk 



2(E Bfi ± % 



IB ' 



exp 



2m- 



L 



T osc 
L B,jk 



exp 



r coh 
L B,jk 



(42) 



In Eqs. (l4"Tj) and (j4~2*|) the upper (lower) signs correspond to Am 2 k > (Am| fe < 0). 
oscillation and coherence lengths in (1421) are defined in analogy with Eqs. and ©: 



The 



T OSC 



4ttE 



B,0 



Am 2 jk 



AttE 
Am 2 jk 



t coh 
^B,jk 



4£ 2 



7s|Amf fe | 7s|Am| fc | ' 



(43) 



We see that Eq. (|4ip depends not on the individual energies and widths of all external states 
separately, but only on the combinations Eb,o = Eu t B,o—Eu c ,B,o and 7/3 = 7h,b+7h c ,b- In the 
limit of no neutrino oscillations, i.e. when all Am 2 k = or U a j = 8 a j, Eq. (14ip reproduces the 



no-oscillation result ( J24|) obtained in our calculation of the Mossbauer neutrino production 
and detection rates treated as separate processes. 

If the localization condition |Am| fc | <C 2a 2 is satisfied for all j and k, as it is expected to 
be the case in realistic experiments, one can pull the generalized Lamb-Mossbauer factor out 
of the sum in Eq. (14 ip and replace the localization exponentials by unity, which yields 



To B 
AtxL 2 



Y S Y D exp 



a: 



<-.) 



2 \U e k\ 2 Ijk ■ 



(44) 



Here m is an average neutrino mass and Ij k is defined in Eq. (1A4j) . In realistic situations, 
it is often sufficient to consider two-flavour approximations to this expression. Indeed, at 
baselines L ~ 10 m which are suitable to search for oscillations driven by #13, the "solar" 
mass squared difference Am^ is inessential, whereas for longer baselines around L ~ 300 m, 
which could be used to study the oscillations driven by the parameters Am^ and 6*12, the 
subdominant oscillations governed by Am^ and #13 are in the averaging regime, leading to 
an effective 2-flavour oscillation probability. In both cases one therefore needs to evaluate 



/] \Uej\ 2 \Uek\ 2 Ijk — 
j,k=l,2 

c 2 s 2 /Ati 



(75 + 7d)/27t 



(FT W \2 1 (75+7d) 2 
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c 2 s 2 



(E S: o — E Dfi ) 2 + 



(7S+7d) 2 



{A (S) _ A ( D)) ^ Esfi _ Ed o){is _ 1d) + 1 



(7S+7d) 2 " 



E 



5,0 



E D .o + i 



IS -Id 



+ C.C. 

(45) 
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where (B = S,D) denotes the value of A - k corresponding to the appropriate fixed 

Am 2 k = Am 2 (which is defined here to be positive, i.e. Am 2 = |Am 31 | 2 or Am^), s = sin9 
and c = cos 8, with 9 being the relevant two-flavour mixing angle. 

As in the full three-flavour framework, in the absence of oscillations, i.e. for Am 2 = or 9 = 
0, Eqs. (jHJ) and fj4"5l reproduce the no-oscillation rate of Eq. (I24"j) . With oscillations included, 
the first line of Eq. ( l45l) factorizes into the Lorentzian times the v e survival probability, which 
in general contains decoherence factors. Such a factorization does not occur in the second line 
because the first term in the numerator in the square brackets is not proportional to 7s + 7d- 
This term, containing a product of three small differences, is typically small compared to the 
other terms (at least when the Mossbauer resonance condition l-E^o ~ -^d.oI (7s + 7d)/2 is 
satisfied). Still, it is interesting to observe that a naive factorization of T into a no-oscillation 
transition rate and the v e survival probability is not possible when this term is retained. 

In all physically relevant situations, however, the whole second line of Eq. (1451) is negligible 
because so is A^ — A^ D \ Retaining only the contribution of the first line in Eq. ( j4"5l) . from 
Eq. (HI) one finds 



, r B 

' 4vrL 2 



Y S Y, 



s y d exp 



E s,o 



VTir 



(75 + 7d)/2tt 



2s 2 c 2 



J (E, 



5,0 



E 



D,0) 



+ 



L + e" aoL ) cos 



4 

/ Am 2 L 



(46) 



where as,D = {Am 2 / AE 2 )^s,d, so that exp[— ots,DL] = exp[— L/L 2 ^] are the decoherence 
factors (cf. Eqs. ( j42|) and ()43l)). For realistic experiments, one expects the oscillation phase 
(Am 2 /4E)L to be of order unity, so that as^L ~ Js,d/E ~ 1CT 15 , and decoherence effects 
are completely negligible. The second line in Eq. (I46p then yields the standard 2-flavour 
expression for the v e survival probability. 

As we have already pointed out, the contribution of the second line in (145 p to T is of 
order (e~° sL — e~ a ° L ) and therefore completely negligible. It is interesting to ask if there are 
any conceivable situations in which the decoherence exponentials in Eq. ( j46l) should be kept, 
while the contribution of the second line in (T4"5"l) can still be neglected. Direct inspection of 
Eq. ( T4"5l) shows that this is the case when \Es,o — Ed,o\ ^ |7s + 7d| with o>s,dL > 1 and 
lets — oid\L <C 1. 



4.3. Homogeneous line broadening 

Homogeneous line broadening is caused by various electromagnetic relaxation effects, in- 



cluding interactions with fluctuating magnetic fields in the lattice [15J, [16|. Unlike inhomoge- 
neous broadening, it affects equally all the emitters (or absorbers) and therefore cannot be 
taken into account by averaging the unperturbed transition probability over the appropriate 
energy distributions of the participating particles, as we have done in the previous subsection. 
Instead, one has to modify already the expression for the amplitude. Since the homogeneous 
broadening effects are stochastic, a proper averaging procedure, adequate to the broadening 
mechanism, has then to be employed. For the conventional Mossbauer effect with long-lived 
nuclei, a number of models of homogeneous broadening was studied in fl6l . 38 . 
all the considered cases the Lorentzian shape of the emission and absorption lines has been 
obtained. The same models can be used in the case of Mossbauer neutrinos; one therefore 
expects that in most of the cases of homogeneous broadening the overall neutrino production 
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- propagation - detection rate will also have the Lorentzian resonance form, i.e. will essen- 
tially coincide in form with Eq. (141 p . or with its simplified version in which the difference 
between A^} and is neglected. A notable exception, which we consider next, is the 
homogeneous broadening due to the natural linewidth. As we shall see, this case is special 
because the time interval during which the source is produced is small compared with the 
tritium lifetime. 



4.4. Neutrino Mossbauer effect dominated by the natural linewidth 



Although in a Mossbauer neutrino experiment with a tritium source and a 3 He absorber 
inhomogeneous broadening as well as homogeneous line broadening different from the natural 
linewidth are by far dominant, we will now consider also the case in which the emission and 
absorption linewidths are determined by the decay widths of the unstable nuclei. Even though 
it is not clear if such a situation can be realized experimentally, it is still very interesting for 
theoretical reasons. 

To take the natural linewidth of tritium into account, we modify our expression for the 
amplitude, Eq. (1261) . by including exponential decay factors in the 3 H wave functions. For 
the tritium in the source, this factor has the form exp(— yt/2), describing a decay starting 
at t — 0, the time at which the experiment starts. 4 For the tritium which is produced in 
the detector, the decay factor is exp(— j(T — t 2 )/2), where ti is the production time and T 
is the time at which the number of produced 3 H atoms is counted. Note that 7 here is the 
total decay width of tritium, not the partial width for bound state beta decay. Since we are 
taking into account the finite lifetime of tritium, we also have to restrict the domain of all 
time integrations in A to the interval [0,T] instead of (—00, 00). We thus have to compute 
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16 



(1 +7 5 )7^e,D 



(47) 



with the same notation as in Sec. 14.11 This form for A can also be derived in a more rigorous 
way using the Weisskopf-Wigner approximation 41, |42|, 31, 43], as shown in appendix O 



4 This is also supposed to be the time at which the number of 3 H atoms in the source is known. It is 
assumed that the source is created in a time interval that is short compared to the tritium mean lifetime 
7 _1 = 17.81 years. 



16 



After a calculation similar to the one described in Sec. 14.11 we find for the total probability 
for finding a tritium atom at the lattice site in the detector after a time T: 
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In the derivation, which is described in more detail in appendix [HI we have neglected the 
energy dependence of the generalized Lamb-Mossbauer factor and of the spinorial terms, 
approximating them by their values at E = ^(Es,o + -Ex),o)- Furthermore, we have expanded 
the oscillation phase around this average energy. These approximations are justified by the 
observation that these quantities are almost constant over the resonance region. 



In Eq. (SHD the quantity Vj = (E< 



-m 



2\l/2 



/ E denotes the group velocity of the jth neutrino 



mass eigenstate, and the generalized Lamb-Mossbauer factor is parameterized in the by now 



familiar form with p"^ m = E — max(m 



m 



Moreover, we have defined the quantity 



L L 
T jk = mm T , T 



(49) 



which corresponds to the total running time of the experiment, minus the time of flight of 
the heavier of the two mass eigenstates Uj and v^. The appearance of the step-function factor 
9 (Tjk) in Eq. (l4*Hl) is related to the finite neutrino time of flight between the source and the 
detector and to the fact that the interference between the jth and fcth mass components 
leading to oscillations is only possible if both have already arrived at the detector. As in 
Sec. 14.21 decoherence exponentials appear, containing the characteristic coherence lengths 



T coll 
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(50) 



In the approximation of ultra-relativistic (or nearly mass-degenerate) neutrinos, this becomes 

AE 2 



j coh 



7 1 Am 



(51) 



jk I 



and is thus analogous to Eqs. (jSJ) and ( j43lh 

While the first two lines of Eq. (l4"8j) contain the standard oscillation terms, the generalized 
Lamb-Mossbauer factor and some numerical factors, the expression in the third line is unique 
to Mossbauer neutrinos in the regime of natural linewidth dominance. To interpret this part 
of the probability, it is helpful to consider the approximation of massless neutrinos, which 
implies Vj = 1 for all j and thus U^ k = oo. If we neglect the time of flight L/i)j compared to 
the total running time of the experiment T, we find that the probability is proportional to 



- 7 T 



sin 2 [(E Sj0 - E D)0 ) 
(E Sfi - E Dfi ) 2 



(52) 
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The factor exp(— yt) accounts for the depletion of 3 H in the source and for the decay of the 
produced 3 H in the detector. 

It is easy to see that for 7 = and T — > 00, Eq. ( 13*71) is recovered, except for the omitted 
averaging over the energies of the initial and final state nuclei. In particular, we see that 
in this limit, due to the emerging 5-function, the Mossbauer effect can only occur if the 
resonance energies Es,o and Ed,q match exactly. For finite T, in contrast, the matching 
need not be exact because of the time-energy uncertainty relation, which permits a certain 
detuning, as long as l-E^o - ^r>,o| ^ 1/T. In the case of strong inequality \Esfi — E Dfi \ <C 1/T, 
Eq. fl""""?]) can be approximated by 

T 2 e" 7T /4. (53) 

It is crucial to note that the allowed detuning of E s $ and E D does not depend on 7, 
contrary to what one might expect. Instead, the Mossbauer resonance condition requires 
this detuning to be small compared to the reciprocal of the overall observation time T. 
Therefore, the natural linewidth is not a fundamental limitation to the energy resolution 
of a Mossbauer neutrino experiment. There is a well-known analogue to this in quantum 
optics (IH, called subnatural spectroscopy. Consider an experiment, in which an atom is 
instantaneously excited from its ground state into an unstable state \b) by a strong laser 
pulse at t = 0. Moreover, the atom is continuously exposed to electromagnetic radiation 
with a photon energy E, which can eventually excite it further into another unstable state 
\a). If, after a time r, the number of atoms in state \a) is measured, it turns out that 
the result is proportional to 1/[(E — AE) 2 + (y a — 7 b ) 2 /4] rather than to naively expected 
1/[{E — AE) 2 + (7 a + 7fe) 2 /4], where AE is the energy difference between the two states, and 
7 a , 7b are their respective widths. In our case, the state \b) corresponds to a 3 H atom in the 
source and a 3 He atom in the detector, while \a) corresponds to a 3 He atom in the source 
and a 3 H atom in the detector. The initial excitation of state \b) corresponds to producing 
the tritium source and starting the Mossbauer neutrino experiment, and the transition from 
I b) to I a) corresponds to the production, propagation and absorption of a neutrino. Since 
the difference of decay widths 7 a — 7^ vanishes for Mossbauer neutrinos, 5 we see that 7 does 
not have any impact on the achievable energy resolution, in accordance with Eq. (15*21) . Note 
that this is only true because the source is produced at one specific point in time, namely 
t = (more generally, during a time interval that is short compared to the tritium lifetime). 
In a hypothetical experiment, in which tritium is continuously replenished in the source, an 
additional integration of V over the production time would be required, and this would yield 
proportionality to l/[(Es,o — Ep^) 2 + j 2 ], in full analogy with the corresponding result in 
quantum optics jlij . 

The T-dependence of V, as given by Eq. ( 1531) can be understood already from a classical 
argument. If we denote the number of 3 H atoms in the source by N$ and the corresponding 
number in the detector by No, the latter obeys the following differential equation: 

N D = -NsNoPee^X ~ jN D ■ (54) 

Here P ee is the v e survival probability, iVo is the number of 3 He atoms in the detector, which we 
consider constant (this is justified if the number of 3 H atoms produced in the detector is small 



5 We assume that the tritium nuclei in the source and detector have the same mean lifetime. 



18 



compared to the initial number of 3 He), and cr(T) is the absorption cross section. It depends 
on T because, due to the Heisenberg principle, the accuracy to which the resonance condition 
has to be fulfilled is given by T" 1 . If we describe this limitation by assuming the emission 
and absorption lines to be Lorentzians of width 1/T, we find that for j-E^o — Ed,q\ <C T _1 
the overlap integral is proportional to T, so that we can write a = sqT with sq a constant. 
Using furthermore the fact that Ns = Ns,o exp(— 7T), the solution of Eq. (1511) is found as 

Nd= 8^ Te • (55) 
This expression has precisely the T-dependence given by Eq. 



5. DISCUSSION 

Let us now summarize our results. We have studied the properties of recoillessly emitted 
and absorbed neutrinos (Mossbauer neutrinos) in a plane wave treatment (Sec. |2|, in a 
quantum mechanical wave packet approach (Sec. [3]) and in a full quantum field theoretical 
calculation (Sec. HJ). The plane wave treatment corresponds to the standard derivation based 
on the same energy approximation. We have pointed out, that for Mossbauer neutrinos this 
approximation is justifiable, even though for conventional neutrino sources it is generally 
considered to be inconsistent. The wave packet approach is an extension of the plane wave 
treatment, which takes into account the small but non-zero energy and momentum spread of 
the neutrino. Finally, the QFT calculation is superior to the other two, in particular, because 
no prior assumptions about the energies and momenta of the intermediate-state neutrinos 
have to be made. These properties are automatically determined from the wave functions of 
the external particles in the source and in the detector. For these wave functions we used 
well established approximations that are known to be good in the theory of the standard 
Mossbauer effect. 

In all three approaches that have been discussed, we have consistently arrived at the pre- 
diction that Mossbauer neutrinos will oscillate, in spite of their very small energy uncertainty. 
The plane wave result, Eq. (17j) , is actually the standard textbook expression for the v e sur- 
vival probability, and Eqs. (jHJ), (141j) and (1481) are extensions of this expression, containing, 
in particular, decoherence and localization factors. We have found that these factors cannot 
suppress oscillations under realistic experimental conditions, but are very interesting from 
the theoretical point of view. 

Let us now compare the results of different approaches. First, we observe that the de- 
coherence exponents in our QFT calculations are linear in L/L coh , while in the quantum 
mechanical result, Eq. (jSJ), the dependence is quadratic. This behaviour can be traced back 
to the fact that Gaussian neutrino wave packets have been assumed in the quantum mechan- 
ical computation, while in our QFT approach we have employed the Lorentzian line shapes, 
which are more appropriate for describing Mossbauer neutrinos. The linear dependence of the 
decoherence exponents on L/L coh in the case of the Lorentzian neutrino energy distribution 



has been previously pointed out in [31 



Even more striking than the differing forms of the decoherence exponentials is the fact 
that a localization factor of the form exp[— | Am^,|/2<jp] is present in Eqs. ( HT1) and (1151) . while 
the localization exponentials disappear from Eq. flSJ in the limit £ — > which is relevant for 
Mossbauer neutrinos. This shows that the naive quantum mechanical wave packet approach 
does not capture all features of Mossbauer neutrinos. In particular, it neglects the differences 
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between the emission (and absorption) probabilities of different mass eigenstates, which ef- 
fectively may lead to a suppression of neutrino mixing. In realistic experiments, however, 
this effect should be negligible. 

Another interesting feature of the exponential factors implementing the localization con- 
dition in our QFT calculations is that the corresponding exponents are linear in |Am^|, 



whereas the dependence is quadratic (for £ ^ 0) in the quantum-mechanical expression 
This can be attributed to the fact that we consider the parent and daughter nuclei in the 
source and detector to be in bound states with zero mean momentum (but non-zero rms 
momentum). This is also the reason why the <r p -dependence of the localization exponents in 
Eqs. (14 ip and 048|) is different from that in the quantum mechanical approach (namely, they 
depend on |Am| fc |/2<Tp rather than | Am'j k \/2pa p ): for the considered bound states, the rms 
momentum is p ~ a p , so that pa p ~ a p . 

One more point to notice is that while the same quantity, the momentum uncertainty a p , 
enters into the decoherence and localization factors in the quantum mechanical formula (JBJ), 
this is not the case in the QFT approach, where the localization factors depend on a p , whereas 
the decoherence exponentials are determined by the (much smaller) energy uncertainty. In 
the case of natural line broadening this energy uncertainty is given by the 3 H decay width 7, 
while in all the other cases it is given by the widths of the neutrino emission and absorption 
lines, which are determined by the homogeneous and inhomogeneous line broadening effects 
taking place in the source and detector. 

The QFT results of Eqs. (14ip and ( 1481) describe not only the oscillation physics, but 
also the production and detection processes. These results can thus also be used for an 
approximate prediction of the total event rate expected in a Mossbauer neutrino experiment. 
Both expressions contain the Lamb- Mossbauer factor (or recoil- free fraction), which describes 
the relative probability of recoilless decay and absorption of neutrinos. Moreover, they contain 
factors that suppress the overall process rate T unless the emission and absorption lines 
overlap sufficiently well. In the case of inhomogeneous line broadening (Sec. 14.21) as well as 
for homogeneous broadening different from the natural linewidth effect, this is a Lorentzian 
factor, the same as in the no-oscillation rate (1241) . It suppresses the transition rate if the 
peak energies of the emission and absorption lines differ by more than the combined linewidth 
1s+1d- We have, however, found that the factorization of the total rate into the no-oscillation 
rate including the overlap factor and the oscillation probability is only approximate. For the 
hypothetical case of an experiment in which the neutrino energy uncertainty is dominated by 
the natural linewidth 7 (Sec. I4.4j) . we have found that the overlap condition does not depend 
on 7, but is rather determined by the reciprocal of the overall duration of the experiment T. 
Although this result may seem counterintuitive at first sight, it has a well-known analogy in 
quantum optics 44[ and is related to the fact that the initial unstable particles in the source 
are produced in a time interval much shorter than their lifetime. 

Notice that the overlap factors contained in our QFT-based results for the neutrino Moss- 
bauer effect governed by the natural linewidth and by other line broadening mechanisms, 
Eqs. (jUj) and ( 1441) . are two well-known limiting representations of the 5-function, which yield 
the energy-conserving 5-function 5(Es,o — Ed,q) in the limits T — > 00 or 75 + 70 — * 0, re- 
spectively. 6 One can see that in these limits both expressions reproduce, if one sets the v e 
survival probability P ee to unity, the no-oscillation rate (1221) obtained in the infinitely sharp 
neutrino line limit by treating the Mossbauer neutrino production and detection as separate 



6 Eq. (Tl8l) yields T6(Es,o ~ En,o) because it describes a probability rather than a rate. 
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processes. Our QFT results thus generalize the results of the standard calculations and allow 
a more accurate and consistent treatment of both the production - detection rate and the 
oscillation probability of Mossbauer neutrinos. 

To conclude, we have performed a quantum field theoretic calculation of the combined 
rate of the emission, propagation and detection of Mossbauer neutrinos for the cases of 
inhomogeneous and homogeneous neutrino line broadening. In both cases we found that the 
decoherence and localization damping factors present in the combined rate will not play any 
role in realistic experimental settings and therefore will not prevent Mossbauer neutrinos 
from oscillating. 
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Appendix A: DERIVATION OF THE TRANSITION RATE FOR 
INHOMOGENEOUS LINE BROADENING 

In this appendix we give details of the steps leading from Eq. f[3"T|) to Eq. (]41|) when the 
energy densities of the external states pa,b are chosen to have the Lorentzian form, Eq. ( j40l) . 

First, we notice that the Lorentzians are sharply peaked, with their widths much smaller 
than the peak energies, and therefore the energy integrals in ( 137|) get their main contributions 
from the narrow intervals around these peaks. Since the peak energies, as well as their 
differences E^s,o — -Eue.s.o an d En S) o — -Efic.s.o that determine the neutrino energies, are much 
larger than the neutrino masses, one can employ the expansion \J E$ — mf ~ E$ — mf/2Es 
in the exponents. 

Next, we make use of the identity 

dE n dEh -^- J - 5 "^—L f(E„ — Eh) 

k (E^E a , t ,r- + ^(E t -E K0 Y+-i n 



that holds for any function f(E) for which the integrals in ( 1A1I) exist. To apply this for- 
mula to Eq. ( 1371) . we have to extend the domain of the energy integrals from the physical 
region [maxj (m^ ) , oo) to the whole real axis, (—00,00). This is again possible because the 
Lorentzians Pa,b{Ea,b) have very narrow widths and therefore ensure that the unphysical 
contributions are strongly suppressed, the error introduced by the extension of the integra- 
tion interval being of order js(D)/ Es(d),o ~ 10^ 15 . We can thus use Eq. (I A 1 [) to perform two 
of the four energy integrations in Eq. (13 7\i . Of the remaining two, one is trivial due to the 
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Im E 



Re E 



E s , - t¥ 



\ 



E D , -i^f- / 



Figure 2: Integration contours in the complex E plane. 



factor S(Es — -Ex>), so that the expression for T becomes 



47TL2 » (E - E Sfi y + 7l/4 (E - E D ^ + 7 f,/4 



2 |£7 efc | 2 exp 



2i? 2 — m 2 - — ml 1 



(A2) 



Next, we pull the generalized Lamb-Mossbauer factor out of the integral, replacing it 
by its value at E = (Es,q + Ed$)/2. This is justified by the observation that 7s, 7d ~ 
1CT 11 eV <C a p ~ 10 keV, so that the generalized Lamb-Mossbauer factor is nearly constant 
over the region where the integrand is sizeable. We are thus left with the task to compute 
the expression 



ijk 



' dE ^ ^ e"^* (A3) 

(E - £ s , ) 2 + 7 §/4 (E - E Dfi f + 7 2 /4 ' [K6) 



which can be done by integration in the complex plane. The integrand of (1A3j) has four poles, 
two above the real axis and two below, and an essential singularity at E — (see Fig. [2]). To 
circumvent the essential singularity, we choose the integration contour to consist of the real 
axis with a small interval [— e, e] cut out, supplemented by a half-circle of radius e around 
the point E = and closed by a half-circle of large radius. The contribution of the small 
half-circle vanishes when its radius goes to zero provided that we avoid the point E = 
from above when Am 2 fc > and from below when Am 2 fc < 0. Thus, we close the integration 
contour in the upper half-plane for Am 2 fc > and in the lower half-plane for Am 2 fc < 0. The 
contribution from the large half-circle vanishes when its radius tends to infinity because the 
product of two Lorentzians goes to zero as \E\~ A for \E\ — > oo, while the exponential becomes 
unity in this limit. Application of the residue theorem yields now 
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ijk 



2nE 



5,0 



E D ,o ± i 



; is -Id 



(5) 



E 



5,0 



Sn ±i2a±2a 



75^ } 



(A4) 



with the notation from Sec. 14.21 Here the upper (lower) signs correspond to Am^ > 
(Ami < 0). Inserting this result into Eq. ( 1A2I) . we obtain Eq. 



Appendix B: DERIVATION OF THE TOTAL TRANSITION PROBABILITY FOR 
THE CASE OF NATURAL LINEWIDTH DOMINANCE 



In this appendix we describe the derivation leading from Eq. (1471) to Eq. (1481) . The spatial 
integrals in Eq. (1471) are the same as those encountered in Sec. 14.11 and yield the factor 



exp[— p /2a 2 ] exp[ipL]. The time integrals can also be evaluated straightforwardly; however, 



unlike the corresponding integrals in Eq. (1261) . they do not give exact, but only approximate 
energy conserving factors for the production and detection processes. This behaviour can 
be ascribed to the non-zero width of the tritium states and the finite measurement time T. 
To evaluate the three-momentum integral over p, we again employ the Grimus-Stockinger 
theorem and thus find 



-IT 



J 

-i(E s -po)T- 



_ I e i(E D -p )T+%T _ j 



Po ~ E £ 



Po ~ E D 



exp 



2 2 

Po ~ 
2a 2 



where the 4-vector pj is defined as pj = (po, (pi — m 2 ) 1 ' 2 'L/L). The exponential depending 
on a 2 , which will eventually lead to the generalized Lamb-Mossbauer factor, can be approxi- 
mated by its value at E = (Eg + Ed)/2 because 7 <C o v ensures that it is almost constant in 
the region from which the main contribution to the integral comes, namely the region where 
|po — Es\ < 7 and \po — Ed\ < 7- The fact that this region is very narrow also allows us to 
pull the spinorial factors out of the integral and to expand the oscillation phase around E: 



m 2 L 



E 2 - m 2 L 



i—(Po 



E) 



(B2) 



where Vj = (E 2 



2U/2 



/ E. The integral over p can then be evaluated by complex contour 



integration. The denominator has poles at p = E s — i'j/2 and p = E D — 27/2, and the 
relevant terms in the numerator are 



e -i( Es - Po )T-l 7 T _ A f e i(E D -p )T+%TT _ j 



i(P0-E)- 



ipo 

e "i e 



i(E s -E D )T-E± ipo(T+±) -iE s T-iE±-\ 



u j — e 
-v ' s - 

(A) 

L \ ; TP^'T ATP L 1 1 , 



j e 



JT 



(B) 

-ipo{T-fr) iE D T-iE±- + ^T i±p - 

— e j e ■? + e 3 e 



iE± 



(C) 



(D) 



(B3) 



To close the integration contour, we add to the real axis a half-circle of infinite radius. For 
the terms labeled (A), (B) and (D), this half-circle has to lie in the upper half-plane, while 
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for (C) it has to lie in the upper half-plane for T < L/vj, and in the lower half-plane for 
T > L/vj. As the integrand is holomorphic for Im(p ) > 0, only in this last case the integral 
can be non-zero. The residue theorem then yields 



iA =^lJ2 6 ( T - L M) MSMZM ■ exp 



• luU e)D e 



E" 1 



u e ,sl^-l 5 )(p 1 +m j )(l+ 1 5 ) 



.y i /E 2 -m? j L^_ 



-h(T- 



^ e -±(E s -E D )T) 



Ea-E 



D 



UE s -E D ){T-±) 



.(E S -E D )(T-^) 



(B4) 



where now pj = (E,(E* 



m 



2\l/2 



L/L), and 9{x) is the Heaviside step function. The total 



probability for finding a tritium atom at the lattice site in the detector after a time T is 

V = W, (B5) 

where the bar indicates the average over initial spins and the sum over final spins. Apart 
from these spin sums, no integration over the energy distributions of the initial and final state 
nuclei is necessary as long as only natural line broadening is taken into account, because we 
are dealing with transitions between discrete energy eigenstates. A straightforward evaluation 
of Eq. (|B5]) yields Eq. 



Appendix C: WEISSKOPF-WIGNER APPROACH TO THE EFFECTS OF THE 

NATURAL LINE WIDTH 

In this appendix we use the Weisskopf-Wigner approach (41, 42, 43, 31] to derive Eq. 047p . 
which has been the starting point for our discussion of Mossbauer neutrinos in the regime 
of natural linewidth dominance. In particular, our aim is to substantiate the arguments 
dictating the form of the exponential decay factors by an explicit derivation. 

We can write the Hamiltonian of the system as H = Hq + e lHot Hie~ lHot , where Hi is the 
interaction-representation weak interaction Hamiltonian and Ho is the remainder. In general, 
we will not treat Hi as a perturbation since we are ultimately interested in the depletion of 
unstable states over time, which cannot be adequately described in a perturbative approach. 

One can write an arbitrary state as \ip{t)) = ^2jCj(t)\(j>j), where \<j)j) are the eigenstates 
of Hq. The Schrodinger equation then gives the evolution equations for the coefficients Cj(t): 

ic j {t) = Y,(<l>3\Hi\<Pk)c k {t). (CI) 

k 

For our purposes it will be convenient to slightly modify the notation and classify the different 
states according to their particle content, as shown in Table HI 3 H and 3 He denote the two 
types of atoms in the experiment, and the index S or D shows whether the respective particle 
is initially localized at the source or at the detector. For those states for which we have 
written the electron participating in the reaction and the 3 He + ions separately, we imply 
that the electron may be either free or in an atomic bound state, while for the other states 
only bound electrons are considered. The upper index (i) stands for the initial state, the 
indices (1) through (6) denote intermediate states, and (/) stands for the final state, after 
the decay of the source particle, the absorption of the emitted neutrino in the detector and 
the decay of the produced tritium. The lower indices stand for the various quantum numbers 
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Particles 




Energy 


Coefficient 


State vector 


Initial state 




3 He D 


E® 






Intermediate states 


3 He£, u s , eg 
3 H 5 

3 He£, u s , eg 

3 He s 

3 H 5 

3 He£, u s , eg 


3 He D 
3 H D , u D 
3 ~H-D, v D 
3 H D 

3 HeJ, u D e D , u D 
3 HeJ, u D e~ D , u D 


ft, 

E^ 

h kl 
p(6) 


f! 

ft, 

c (4) 

c (5) 
c (6) 

jkl 


l^ 2) ) 

l r ft, ' 

i$> 

1*1?) 

i*SS> 


Final state 


3 He s 


3 HeJ, u D e D 


4 n 


^ 





Table I: Classification of the states appearing in a Mossbauer neutrino experiment, 
of the particles; for example, j encodes the momenta and the spins of z/$ and eg, and the 



information whether e s is bound or free. 



The evolution of the system is governed by the interaction Hamiltonian Hi = Hg + H D 



Hp + H.c, where 



H+ = Jd 3 x^=G F cose c ( 3 Ee\J^E)^sl,(l-l 5 )^, (C2) 
H D = J d 3 x^=G F cose c ( 3 R\J^\ 3 Be) ^(1-^)^0, (C3) 
H+ = yrf 3 x-J=G F cos^ c ( 3 He|^| 3 H)^ eiS7M (l-7 5 )^. (C4) 

The Hermitian conjugates of these operators will be denoted Hg , Hp and Hp. The Hamil- 
tonians Hg and Hp describe tritium decay in the source and detector respectively, whereas 
Hp describes the v e capture in the detector. Although the Hamiltonians are essentially re- 
lated by Hg = Hp = Hp, we will treat them as distinct operators throughout this appendix 
to keep our derivation more transparent and more general. For the matrix elements of the 
transitions, the following relations hold: 



<0«>li/ s -l 


= ( 


4* 






(<P (<) \H + D \ 


4 2) > = < 








(4> m \fl- D \ 


tf>) = < 


4» 







(C5) 



and similarly, 

B« _ e {1) - E m - E {3) - E {5) - E {6) 

j ~ k ^jk — ^kl ^jkV 

£;(•) - B f = E f - E®, (C6) 
They follow from the fact that the corresponding processes differ only by the spectator 
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particles. The evolution equations for the system are 

- E^i^i^) 4 1} + E<<^i<4 2) > 4 2) , (C7) 



2C V 



A- 



i<f = < W|*">) c"» + £>?'l ff JI$'> cf, (C8) 

4 2) = wPiflBi^) c» + B^'i^f) 4? + E<*S*i^i*2 , > 4?, (c») 

B* (4) i^ W w > 4 1} + £^ (4) i#M 2) > 4 2) + E^ (4) i^ (/) > 4'\ ( cn ) 

j k I 



lc jk ~ Wjk 



MW -<<WI^>4? + <ijSl^§>^, (C13) 

E^ (/) i^M?> 4? + (0! /} i^i0 (4) ) c (4) . (ci4) 



We treat all processes that occur within the source or within the detector non-perturbatively, 
while first-order perturbation theory will be used for processes that require the propagation of 
a neutrino between the source and the detector. This second kind of transitions is suppressed 
due to the smallness of the solid angle at which the detector is seen from the source. Con- 
sequently, we include only the respective forward reactions (i.e. those proceeding downward 
in the scheme of Tabled]), but neglect the feedback terms, which would otherwise appear 
in the equations for Cj , cjj. , and cjy . The feedback of \4>\^) to \4>^) is included because 
the production of both states from the initial state requires a single neutrino propagation 
between the source and the detector. The sums in Eqs. (IC7I) - (IC14p symbolically denote the 
summation over the relevant discrete indices and integration over the continuous variables. 

The initial conditions for the equation system ( 1C7I) - ( 1C14I) are given by c w (0) = 1 with 
all other coefficients vanishing at t = 0. 

Our ultimate goal is to solve the evolution equations for c^'(t), which determines the 3 H 
abundance in the detector at time t. It is convenient to first consider the closed subsystem 
formed by Eqs. §C7§, Q, (jCSj) . (ICTO . (ICT21 and 05ft . which we solve from the bottom 
upwards. We start by integrating Eq. flC13j) to obtain an expression for c^, which we then 
insert into Eq. (]C12p . This yields 

4?(*) = <4?i^(*)l4 2) )4 2) (*) 



— % 

3 



E (<f>$\ H s(t)\<f>%) i&HUti)^) c$(h) . (C15) 

„• JO 



Consider first the second term, which describes the effect on |<^jy ) of its decay into 
Following the Weisskopf-Wigner procedure as described in [43[ , we split the quantum numbers 
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indexed by j into the energy E^ and the remaining parameters (3. Denoting the density 
of states (the number of states per unit energy interval) by p(E<®,/3), one can make the 
replacements 

l*$> - \<$\ E^\P) , E ^ E / dE^p(E^\p) (C16) 

3 P 

in the second term of Eq. (1C15P , which gives 



-i JdE®K(E®) jy^e^-^-^c^ih). (C17) 

Here we have explicitly written down the time dependence of the matrix elements and intro- 
duced the quantity 

K(E®) = J2 \(&Hs(0)\^;E^,^)\\(E^,P), (C18) 

which is a smooth (non-oscillating) function of energy. More specifically, K(E^) represents 
a broad bump of width 0(mw), so that a non-negligible contribution to the energy integral 
in (1C17I) can only arise if t — 1\ < 1/mw- Otherwise, the integrand is fast oscillating and the 
integral is strongly suppressed. Therefore, we can to a very good accuracy replace c^(ii) 
by c^(t) in Eq. ( 1C17I) and pull it out of the integral over t\ (we assume that cf t (t) is 
approximately constant over time intervals of order This assumption will be justified 

a posteriori by inspecting the obtained expression for cjy (t)). For t ^> 1/mw we thus obtain 



J dE® K{E®) f dh e^ L£(6) »^ 



-ic 



(5), 
kl 



(t) dE^ 



^(E^-E^+ipf—^ ) 



2 



-M?+^)4?(*)> (ci9) 



where 



7 = 2kK(E®) , 5E = P ( dE™ % {E _ J , (C20) 

J E kl E^ 

and P denotes the principal value. As follows from the definition of the function K(E) in 
Eq. ( 1C18I) and Fermi's golden rule, 7 is just the decay width of 3 H in the source. The quantity 
5E is the mass renormalization of the particles forming |</>^). From now on, we will omit 5E 
and similar quantities in subsequent formulas, assuming that they are already included in the 
definition of the physical masses of the involved particles. The formal solution to Eq. ( 1C15I) 
is 

= ~i fdh (^iHUm^e-^-^cfih) (C21) 
Jo 

+ H) 2 £ fdt x P dutf$\Hg{t x )\4$ <<|fl£(f a )l$> z- Ht - tl] ■ 



Jo 
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By a similar argument, we obtain from Eq. (jClOP : 



- i fdti {4>tH\H^{t)\<t>%) (<P%m(ti)\<P { S) 4?(*0 , (C22) 



' (3) I 



;(!)\ J 1 ) 



(3) | 



where the decay width of 3 H in the detector, 7, has been defined in analogy with Eq. (1C20I) . 
We will now show that the last term of Eq. (1C22I) can be neglected. To this end, we insert 
in it the expression for cj^ (t) from ( 1C21I) . which yields 



-i) 



dt, f ^K^I^Co)!^?)! 2 ^ 6 ^^ 



> 



n(2) U5) 



o 3 E /* 2 ^3 i<0£ 2) i^(o)|^>i 2 (^Sl^c*!)!^) ^<^£1 ) l^fe)]^-% i > 

, Jo Jo Jo 



A E ^>^ cf k {t 3 ) 



(C23) 



Here we have used Eqs. (105|) and (IC6p . We will show now that the first term of Eq. (IC231) 
can be neglected; a similar argument can be used to justify the neglect of the second term. 
Replacing the index I by E^ and (3 in analogy with Eq. flC16p . we obtain 



dE^K{E {5) ) / dh 



dt- 



(C24) 



with k(E®) defined analogously to K(E^). As in Eq. flCT7L the energy integral is non- 
negligible only if t — t 2 < 1/mw- We see immediately that here this condition also implies 
t — ti < l/m w . Consequently, we may pull out of the integral those terms which remain 
approximately constant over time intervals 0(1/ m\y), which gives 



fi^mm^it) ldE®K{E®) I dt, I dt 2 e 



m w 



7TC 



5(E { ^ -E^) +iP 



E 



(2) 



£(5) 



O 



(-) 

\m w J 



(C25) 



which is negligible compared to the other terms contributing to (t) (cf. Eq. (1C22I) ). This 
result already suggests the general rule that the only transitions which may contribute size- 
ably to the evolution equations are those corresponding to the direct production of the states 
(i.e. production with a minimum number of intermediate steps), and those corresponding to 

direct feedback from a daughter state into its immediate parent state, e.g. from \4>f^} into 
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\4>fk)- All terms corresponding to more complicated interaction chains are negligible. One 

,0). 

jk ■ 



can now solve Eq. (1C23[) for " 



Jo 

- i f dt^Htih)^) e-^- 4l) cf\h) . (C26) 
Jo 

Next, we plug our expressions (jC26D and flC2H for cf k ] and cJJ into Eq. (IC9l : 

i4 2> (*) = (^VdWI^V^) -^4 2) (t) -t|cg?(f) 

- i £ (^I^WI^) (^I^OI^) e"^* ft) • (C27) 

We have omitted a term containing the product of (c/>j^ \Hq\4>$), (^fil^sl^fw an( ^ 
(<j>fki\Hp\(j)fk) and thus describing the transition chain \<pf k ) — > — > ) — > |</4^)> 
because this term can be shown to be O^/rriw) by an argument similar to the one we used 
for Eq. (IC251 The formal solution to Eq. IjflgTj) is 



(t) = -i / cZ*! (4 2) l^(ti)|<^ (i) > e -^*-*i)-i^*-*i) c«(t!) (C28) 
Jo 

+ H) 2 fdh f X dt 2 (4 2) |^(^)l^> (^l^foOl^) e-H-*K^) ef (t 2 ) . 
j </o Jo 

We now proceed to Eq. ( 1C8I) : 

<^ 1) (t) = (^ 1) |fl+(*)|^> C W(*) 

i £ fdh (4>?\H%(t)\4!$) {<t>fi\H- D {t x )\ct>f) e-^(^) <f (t x ) . (C29) 



i 

k 



(2) (2) (3) (1) 

The contributions coming from c k through the transition chain \4> k ) — ► |0k ) — ► |0J ) are 

again omitted as being 0(j/mw)- The term containing {(j>j \H^{t)\<pf k ) (4>^\H^(ti)\<p^) 

describes the direct feedback from \4>f k ) to \4>f ), but since the transition \(f)j ) — * \<f>fk) does 
not occur spontaneously, the corresponding decay width is zero. Indeed, when applying the 
Weisskopf-Wigner procedure, we see that the resulting 5-function under the energy integral 
is zero for all allowed energies. Thus, the second term in Eq. (1C29I) is negligible, and the 
equation is solved by 



#\t) = -i I {4^\H+(h)\^) c®(h) • (C30) 



We can insert this expression, together with c k 2 \t) from Eq. ( 1C28I) . into the equation for 
c^(t), and find 

c «(t)= e -^, (C31) 
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up to a term suppressed by •y/m w . The closed-form expressions for cp (t), (4 (t), c^(t), 
cjy (t), and cf^{t) are then 



4 2) w 



dt! (^ 2) |^(^)|(/»W>e-i 



(0\ „- 31*1 



L JO 



J U0 



(»)\ p -|t(*-*i) 



4?(i) = H) 2 / *i / ^2(Cl^(ti)lC>(Cl^(i2)|^>e-f 
Jo Jo 



d 5 )\ 



d 2 )\ /^( 2 )| 



(0\ -57*-57(tl-t2) 



dt! (^ 1) |F+(t 1 )|^))e-^ 

*i r^(0i?|^(ti)l4 2) ) (0i 2) |^fe)l0 (i) )e-^- t2) 



'0 JO 

In the expressions for c^(t) and c^(t), we have used the identity 

dt\ I dt 2 = / dt 2 / dti- 
Jo Jo Jt 2 



(C32) 
(C33) 
(C34) 
(C35) 



(C36) 



(C37) 



Eqs. (103 lj) - (IC36P show that all coefficients are slowly varying over time intervals of order 
1/mw, which provides the a posteriori justification for pulling them out of the time integrals 
when applying the Weisskopf-Wigner procedure. 

We have now all the ingredients required to solve for c^\t). We insert Eqs. ( 1C32I) . ( 1C33I) 
and ( 1C14I) into Eq. (ICllh . neglect the Ofi/mw) contribution from the reaction chain \4>fi) — > 



df) 



) — > \4>^), and apply the completeness relations 



(C38) 



to dispose of the sums over j and k and of the intermediate bra- and ket-vectors in the 
products of matrix elements. This leads us to the main result of this appendix, 



c^\t) = {-if / dh 

Jo Jo 



dh (4> {4) \ 



(C39) 



We see that (t) is given by the time-ordered product of the two interaction Hamiltonians, 
supplemented by the classically expected exponential decay factors. After inserting the ap- 
propriate expressions for Hg and H^, finally setting 7 = 7 and applying the Feynman rules, 
Eq. (rail leads directly to Eq. gTD of Sec. [OJ 
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For completeness, we also give the expression for c\ f \t) 



,(/) 



it) 



-z 



dti 



dU 



t2 



dU 



<.(/)i 



^+(t 1 )^(t 2 )e-^- t2 )F+(t 3 )e-i 



7*3 



+ HUti)H+(t 2 )e- l ^Hv(t 3 )e 



+ H+{t 1 )e~*'*H+(t a )HB{t a )e-* 



|0«>. (C40) 



Note that an alternative way of solving Eqs. (IC7I) - fl014j) is to exploit the fact that, in 
the closed system formed by Eqs. tfUT]), (jUS]) . flH9|) . (lUTO]) . (1UT21 . and f!UT3|) . the processes 
in the source and those in the detector can be separated by using a product ansatz for the 
coefficients c. Once this subsystem is solved, and c, can be computed as above. 
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